function Output=TempFun_ConsDistIRF(PP,SS,MODEL,EquJac,IRF)

List_Shock      =   {'Eps_Eta_M_dom','Eps_Eta_M_ext','Eps_Eta_Y_H'};

% IRF of Consumption Distribution
IRF_ConsDist    =   struct();

Cons_Diff       =   struct();
Cons_VarDec     =   struct();
DevUIP          =   struct();
ir_DevUIP       =   struct();
for ii_Irf=1:length(List_Shock)
    SH                  =   List_Shock{ii_Irf};
    IRF.(SH)            =   IRF_Aug(PP,SS,MODEL,IRF.(SH));
    IRF_ConsDist.(SH)   =   IRFAdd_ConsDist(PP,SS,MODEL,EquJac.HH_Aux,length(IRF.(SH).C),IRF.(SH));
    
    C_byFI_dom          =   IRF.(SH).AggStat(1,:);
    C_byFI_ext          =   IRF.(SH).AggStat(2,:);
    C_byRI_N            =   IRF.(SH).AggStat(3,:);
    C_byRI_H            =   IRF.(SH).AggStat(4,:);
    C_byNA_poor         =   IRF.(SH).AggStat(5,:);
    C_byNA_rich         =   IRF.(SH).AggStat(6,:);
    
    [~,MaxInd]          =   max(abs(IRF.(SH).C));
    C_AggPeak           =   abs(IRF.(SH).C(MaxInd));
    
    Cons_Diff.(SH)      =   [C_byRI_N-C_byRI_H; ...
                             C_byFI_dom-C_byFI_ext; ...
                             C_byNA_poor-C_byNA_rich]/abs(C_AggPeak);
    Cons_VarDec.(SH)    =   [IRF_ConsDist.(SH).VarDec.RI(1,MaxInd);...
                             IRF_ConsDist.(SH).VarDec.FI(1,MaxInd);...
                             IRF_ConsDist.(SH).VarDec.NA(1,MaxInd)];
    Cons_VarDec.(SH)    =   100*(1-Cons_VarDec.(SH));
    
    DevUIP.(SH)         =   100*SubFun_MaxAbsVal(IRF.(SH).DevUIP);
    ir_DevUIP.(SH)      =   SubFun_MaxAbsVal(IRF.(SH).ir./IRF.(SH).DevUIP);
end

Output  =   [ir_DevUIP.Eps_Eta_M_dom;...
             ir_DevUIP.Eps_Eta_M_ext;...
             Cons_Diff.Eps_Eta_M_dom(2);...
             Cons_Diff.Eps_Eta_M_ext(2);...
             Cons_VarDec.Eps_Eta_M_dom(2);...
             Cons_VarDec.Eps_Eta_Y_H(1);...
             Cons_VarDec.Eps_Eta_M_ext(2)];